Linking rigidity transitions with enthalpic changes at the glass transition and the 

fragility of glass-forming liquids 
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A low temperature Monte Carlo dynamics of a Keating like oscillator model is used to study 
the relationship between the nature of glasses from the viewpoint of rigidity, and the strong-fragile 
behaviour of glass-forming liquids. The model shows that a Phillips optimal glass formation with 
minimal enthalpic changes is obtained under a cooling/ annealing cycle when the system is optimally 
constrained by the harmonic interactions, i.e. when it is isostatically rigid. For these peculiar 
systems, the computed fragility shows also a minimum, which demonstrates that isostatically rigid 
glasses are strong (Arrhenius-like) glass-forming liquids. Experiments on chalcogenide and oxide 
glass-forming liquids are discussed under this new perspective and confirm the theoretical prediction 
for chalcogenide network glasses. 

PACS numbers: 61.43.Fs-61.20.-x 



The question of the liquid to glass transition and the 
nature of its characteristic temperature T g have received 
a huge amount of interest in the recent years with a 
special emphasis on the dynamic properties of the glass- 
forming liquid [H-Q- Experimentally, the slowing down 
of the dynamics is mostly tracked from viscosity (or 
structural relaxation time r Q ) measurements. The be- 
haviour of these quantities with inverse temperature not 
always displays an Arrhenius-like (or strong) behaviour. 
When properly rescaled with 1/T g in a semi-log plot, the 
viscosity or relaxation time can indeed display a variety 
of different behaviours usually quantified by a fragility in- 
dex M Q which characterizes the steepness of the slope 
of the relaxation time near the glass transition: 
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M ranges typically between 16 for the strong (Arrhcnius 
behaving) silica liquid, whereas M is about 340 for the 
fragile ortho-terphenyl which changes its viscosity by 10 
orders of magnitude over only 50 K temperature change. 

There have been various efforts to connect the liq- 
uid fragility to some easily measurable quantities in the 
glassy state such as compressibility [6] or Poisson ra- 
tio pj, or to the out-off equilibrium behaviour At 
a somewhat more microscopic level, combined effects of 
the structure and the local mechanical behaviour arising 
from the interaction potential has to act on the dynam- 
ical and calorimetric behaviour at the glass transition 
and thus on the fragility. From a theoretical viewpoint, 
one may therefore wonder how the effect of the potential 
strength in the glassy state on the relaxing behaviour at 
higher temperatures can be determined from simple but 
insightful models. 

The present Letter attempts to address this basic is- 
sue by following the Monte Carlo dynamics of a harmonic 



Keating-like oscillator model that mimicts the elastic fea- 
tures (flexible, isostatic, stressed rigid) of the glass. The 
solution of the model shows that isostatic (i.e. optimally 
constrained or intermediate) glasses exhibit a minimum 
in the energy change during a cooling-annealing cycle 
through the glass transition, independently of the applied 
cooling rate, a result that matches exactly the Phillips 
optimal glass condition [s|. Furthermore, the present 
minimum is correlated with the minimum obtained in 
the fragility index. These findings are successfully con- 
fronted with experimental data on chalcogenide network 
glasses. One can thus conclude that isostatic glasses are 
strong glass-forming liquids, a result that opens new per- 
spectives for the description of glass-forming liquids from 
the viewpoint of the mechanical behaviour of the glass. 
And since rigidity can be also tuned in soft solids and 
colloids (lOj |. it may provide a general clue for an im- 
proved understanding of the dynamic properties leading 
to a strong viscosity behaviour. 

The elastic nature of network glasses with changing 
connectivity (or mechanical constraints n c derived from 
Lagrange-Maxwell counting [1 if) can be modelled in the 
framework of rigidity theorv[ia| . [l3| using a Keating 
potential that represents a semi-empirical description of 
covalent bond-stretching (BS) and bond-bending (BB) 
forces. From this description, it appears that the num- 
ber of zero frequency (floppy) modes / behaves as 3 — n c . 
The glass composition at which one has the vanishing of 
/ corresponds to the Maxwell-Lagrange isostatic condi- 
tion n c = 3. We consider a network of N atoms having 
two types of harmonic oscillators with respective density 
/ and n c . The first has a typical frequency w associated 
with the harmonic motion due to the floppy modes, i.e. 
the modes that allow a local distorsion of the network 
with a low cost in energy. The second has a typical fre- 
quency f2 and represents the Keating potential containing 



the constraints imposed by BB and BS forces 1J, [l5| . A 
similar approach has been recently used to describe the 
effect of rigidity on the glass transition temperature . 
The interaction potential V is then given by: 



V = e 



E = \j^fJx* 



+ n c n 2 xf 



(2) 



where u> and Q represent respectively the floppy mode 
frequency and a typical vibrational mode related to BS 
or BB interactions. An inelastic neutron study [17[ of 
the ternary network glass Ge-As-Se provides information 
about the order of magnitude of the floppy mode, and 
the typical BS and BB vibrational frequencies (energies). 
The floppy mode energy is about 4 meV while indenti- 
fied bond-stretching and bond-bending vibrations have a 
respective energy of 31 and 19 meV, i.e. about six times 
more than for the floppy mode energy. Thus one has 

~ 6w and this ratio will be used for the forthcoming 
numerical applications. 

Glass transition: We now build on the approach de- 
veloped by Ritort and co-workers [3] [HI, i.e. an en- 
ergetical move AV is realized on the oscillators. It is 
applied according to the Metropolis algorithm, i.e. ac- 
cepted with probability 1 if the energy decreases, other- 
wise with a probability exp(— f3AV). The position of the 
Keating and the floppy mode oscillators are simultane- 
ously shifted by Xi + Ri/\f~N and Xi + ri/y/N where Ri 
and Ti are random variables with a Gaussian distribution 
having zero mean and respective variance A 2 and S 2 . 

In comparing the nature of both oscillators, one ob- 
viously has iJ > A as motion around the position Xi is 
facilitated by floppy modes, whereas the amplitude Xi 
of the 'Keating' oscillators should be restricted to small 
vibrations, typically a fraction of the interatomic bond 
distance [20(. The transition probability for the change 



in energy AU is : 

p(AF) =/I(nv&-p^ 2 ) 
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Using the Fourier transform representation of the delta 
function, equ. ([3]) can be expressed as a Gaussian in the 
limit N —> oo as: 
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with V a = \ \ju?8 2 + n c Q 2 A 2 } and: 



Q(t) = n c n 2 A 2 < E{t) > +fu 2 8 2 < e(t) > 



(4) 
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FIG. 1: Top panel: Energy V(T) of the system (solution of 
equ. [7} under cooling and annealing (q = ±1 K/s) for three 
selected systems: n c = 2.2 (flexible), n c — 3.0 (isostatically 
rigid), n c — 4.0 (stressed rigid). The two last curves have 
been shifted upwards for a clearer presentation. Broken lines 
represent the equilibrium state T/2. The shaded areas (e.g. 
on n c = 4) serve to quantify the enthalpic changes (H). 



where the mean position of the oscillators have been 
taken as zero and the averages performed over different 
dynamical histories but with same initial condition for 
the ensemble. From the Metropolis rule, one can then 
write the equation for the evolution of the energy: 



dV 
~dt 



xP(x)dx + / xP{x)e 
Jo 
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dx 
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which, applied to the probability distribution (|4|) leads 
to: 



dV 
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with er/c the error function, and 



g(t) = cxp[-(3Vo + p 2 Q(t) er/c 



2pQ{t) - Vo 
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Equation is not closed, because it depends on Q(t) 
through the time dependence (equ. ([5])) of the averages of 
the floppy mode or constraint energies e(t) and E(t). For 
simplicity, we solve equation © in the low temperature- 
long time adiabatic approximation [2l[ where the deriva- 
tive of the energy V iz zero. One then has: 



Q(t) = |(2V ) 



2V V(t) 



(9) 



First, one can solve equation (J7J) to obtain the behaviour 
of the energy V(T). In order to highlight the effect of the 



number of constraints n c on V(T), we work at fixed coo 
ing/heating rate q = T = ±1 K.s~ l and use To = 500 h 
8 = 10A, £1 = 0.3. Figure [1] represents the evolution < 
V(T) for three different values of n c . First, one not( 
that at high temperature, the energy of the system 
equal to T/2 which is also the equilibrium solution < 
equ. {7J, in agreement with the equipartition theoren 
Glassy behaviour (i.e. deviation from the T/2 line) or 
sets at lower temperature, defining a glass transition r< 
gion around 150 K. Annealing from a low temperatui 
end point shows the typical hysteresis behaviour that 
usually manifested in experiment by an enthalpic ove: 
shoot in the heat capacity [22| • From Figure [U one se< 
that the enthalpic change between the cooling and hea' 
ing curves (quantified by the area) depend on the numbc 
of mechanical constraints. Furthermore, these change 
are minimized for a system that is nearly isostaticall 
rigid (n c = 3), i.e. when the trial moves can be on] 
realized on the oscillator with the highest frequency f 
Additional stiffening (i.e. increase of n c ) of the systei 
for n c > 3 leads to a global increase of Vq and thus to a 
increase of the area. 

The area H defined by the difference between the coo 
ing and heating curves can be tracked with the nun 
ber of constraints n c and the corresponding behavioi 
is displayed in Figure [21 It shows that for a fixed cooi- 
ing/heating rate q, certain glass transitions occur with 
minimal enthalpic changes. The model is therefore able 
to reproduce the Phillips phcnomonclogy of ideal glass 
formation when a glass is optimally constraint, i.e. iso- 
static @, [12] • This is experimentally observed on a char- 
acteristic enthalpy AH nr (see also Fig. [5]) extracted from 
complex heat flow measurements at the glass transition 

Dynamics: The study of the dynamics of the system 
can be achieved from the linearization of equation ([7|) in 
the vicinity of the equilibrium value of V(T) that leads 
to a typical relaxation time: 



1 



T =-^nTZV*expV Q /4T 



(10) 



This allows to define a fragility index M from equ. (TT|) 
which is represented in Figure [2] (right axis) as a function 
of the number of mechanical constraints n c . One can 
clearly remark that the fragility behaviour parallels the 
enthalpic changes at the glass transition. A minimum 
is found for both quantities when the number of floppy 
modes / vanishes, i.e. when the glass is isostatic. This 
point corresponds to the location of the flexible to rigid 
transition 12] and the centroid of the intermediate phase 
0- 

Comparison with experiment. Is there any experimen- 
tal correlation between isostatic glasses and strong glass- 
forming liquids ? 

Fairly complete viscosity and calorimetric measure- 
ments on several covalent glass-forming liquids are avail- 
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FIG. 2: Enthalpic changes at the glass transition as a function 
of the number of mechanical constraints n c for different cycles 
at fixed cooling/heating rates q = 0.1 K/s, q — 1 K/S and 
q = 10 K/s. Right axis: fragility of the system. 



able in the literature. We focus on systems that undergo 
a flexible to rigid transition, i.e. oxide and chalcogenidc 
glass in selected composition ranges that lead to a mean 
coordination number of f = 2.4 and n c — 3 [§],[l2j]. Com- 
positional trends in the activation energy for viscosity 
Ea for binary and ternary chalcogenide and oxide liq- 
uids appear in Fig. [3) In the same figure (right axis) are 
also projected the non-reversing enthalpies AH nr of cor- 
responding glasses. The latter quantity provides an ac- 
curate measure of the enthalpic changes that have taken 
place du ring a heating/cooling cycle at the glass transi- 
tion 



22, 23, 25, 26]. One can note that the global min- 
ima in the activation energy Ea coincide at n c ~ 3 with 
those in AH nr for chalcogenide glasses. Together with 
Fig. [2l these data demonstrate the correlation between 
the strong-fragile classification of glass forming liquids 
with the flcxible-intermediate-stressed-rigid classification 
of corresponding glasses. The correlation unequivocally 
shows that Intermediate phase glasses where n c ~ 3 give 
rise to strong liquids, while both flexible (n c < 3) and 
stressed-rigid glasses (n c > 3) give rise to fragile liquids. 

Chalcogenide glasses can be accurately described with 
a harmonic Keating potential [12, 0] ■ Can the present 
conclusions be extended to glass-forming liquids charac- 
terized by potentials different from those shown in equ. 
^ ? The bottom panel of Figure [3] already sketches 
some limitations. In alkali silicate glasses, a large value 
for n c corresponds indeed to the silica-rich compositional 
region, i.e. to systems where the interaction can be fairly 
described by a Keating potential [13]. However, larger 
amounts of alkali ions (i.e. leading to lower n c 's) increases 
the number of more weaker (Coulombic) interactions and 
should cancel the correlation. Weak alkali atom-non- 
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FIG. 3: Comparison between the activation energy for vis- 
cocity Ea and the non-reversing heat flow AH nr at the glass 
transition (right axis) , as a function of the number of mechani- 
cal constraints n c for three selected glass systems: As-Se (data 
from 0], [H|); As-Ge-Se (data from [13|; sodium sili- 
cates (data from [5sJ, [1^]). The vertical broken lines indicate 
the isostatic composition. 



bridging oxygen ionic bonds form, and as T > T g , these 
weaker interactions cease to act as mechanical constraints 
enhancing the alkali-atoms mobilities and contributing to 
the fragility. In the modified oxides as in the H-bonded 



systems [31|, one therefore does not expect the glass- 
liquid correlation to uphold as in the chalcogenides. 

In summary, we have shown that a statistical model 
using a Keating potential with floppy modes was able 
to reproduce the generic features of the glass transition 
from a simple Monte Carlo dynamics. Glasses which are 
optimally constrained (isotatic) give rise to strong glass- 
forming liquids, and are found to display glass transitions 
with few enthalpic changes. Comparison with experi- 
ments shows that the demonstrated relationship holds 
for network glass- forming liquids. Limitations appear for 
systems where the nature of the potential V is obviously 
changed in a more deeper fashion. 

It is a pleasure to acknowledge ongoing discussions 
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